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We have found for the general class of Modified Gravity Models f{R, G) a new instability which 
can arise in vacuum for the scalar modes of the cosmological perturbations if the background is not 
de Sitter. In particular, the short- wavelength modes, if stable, in general have a group velocity which 
depends linearly in fc, the wave number. Therefore these modes will be in general superluminal. We 
have also discussed the condition for which in general these scalar modes will be ghost-like. There is 
a subclass of these models, defined out of properties of the function f{R, Q) and to which the f{R) 
and f{Q) models belong, which however do not have this feature. 



I. INTRODUCTION 



Modified Gravity Models (MGM) have been introduced as an alternative to the quintessence picture for dark energy 
P, 0, H, 3, [^, 0, 0j B @1- The idea is speculative but attractive: maybe the acceleration of the universe is due to a 
modification of the behaviour of gravity at large scales. Gravitation seems to be described by an action which is not 
renormalizable, and the symmetry of the theory, does not restrict enough the possibility for building up an action for 
the g^y variables, the metric tensor components, as already pointed out in 1917 [l3,[iil- However, to build a sensible 
gravitational theory is not an easy task. Some tried to change it to have a renormalizable theory at 1-1oot 



e tneory at i-ioop 12 , 

but in generahthis leads to fourth-order gravity theories and spurious degrees of freedom p^[Ta[l7l.[l8l.[T9l.l20l.El 



13|,[14|, 

22| 



Lovelock [23[ found and studied a special class of scalars, the so called Lovelock scalars, built out of the Riemann 
and metric tensors only. In four dimensions the only such scalars which are not identically zero are three: a constant, 
the Ricci scalar, and the Gauss-Bonnet combination Q. This last scalar is made of quadratic combinations of the 
Riemann tensor, which in turn depends on second derivatives of g^^i,. Naively, one might expect that this term in a 
Lagrangian would automatically lead to four derivatives for g^^. However, all the Lovelock scalars have the property 
(which makes them special among all possible scalar invariants) that the equations of motion still remain of second 
order. This is still true even if these same curvature invariants are coupled to some other matter scalar field. Another 
property of the Gauss-Bonnet term is that, in four dimensions, it can be written as a total derivative [23 |. therefore 
the only way for it to give a contribution to the equations of motion is to couple it to some other field. Indeed, string 
theory seems to predict the existence of such couplings in the low-energy effective action [2^ [2^ [l^, H^, . 

In the aim of building up the action for gravity, since general covariance is not strong enough, we can take a 
bottom-up approach and look for the most general action which can be built out of the Riemann and the metric 
tensor. Therefore recently models have been introduced which use one or both these scalars. These are the so called 

■ H [13, Si , f{Q) i, i m, M M M, m, a 




331, iM, 21, M, M, m 

and f{R,Q) theories [3 





42 




61 




All the Lagrangians written through these 
scalars share the property that they can be rewritten as scalar fields coupled to R and Q respectively [g^ . Therefore 
introducing functions of R and Q adds scalar degrees of freedom only to the theory, eliminating the spurious spin-2 
extra degrees of freedom which are typically ghost-like [H [H, [el, [66| . 

It is well known that f{R) theories can be written as a scalar tensor theory which, through a conformal transfor- 
mation, takes in the Einstein frame the same form of the Lagrangian for a scalar field minimally coupled to gravity 
[stLIgs}. However, the same technique does not prove it useful for the f{ Q) theories. The vacuum structure of these 
theories has been deeply studied for maximally symmetric spaces (20l . |62| . As for the behaviour of the same theories 
on a Friedmann-Lemaitre- Robertson- Walker (FLRW) background some papers appeared [6^, [6^, [6^, but, as we shall 
see later, they only treated a subset of the whole class of /(-R, G) theories. 

In this paper we will give a general study of cosmological perturbations in vacuum for a general J{R, Q) Lagrangian. 
Some papers claimed that they studied the most general MGM [gl, [6^. However, we think that only a subset of 
these theories was considered (among these, the f{R) and the /(G) theories). We present here a new result for the 
most general MGM, that is the dispersion relation for the scalar perturbations changes, and typically scalar modes, 
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if stable, will be superluminal with a speed proportional to the wave number k. This change for the dispersion 
relation is always present except for a special subclass, which actually contains all the models for which cosmological 
perturbation theory was discussed [H, [6^. For such models, we recover the same results found before. 

This new feature is essential to study the vacuum structure for these theories and an additional bound (the speed 
of the short wavelength modes) should be added to the previous one found in the literature of these models. We want 
to stress out that this is a physical property of the general class of these theories, and not a feature of a bad gauge 
choice, or a spurious feature which can be removed by a suitable field redefinition. Although the study of this action 
in the presence of matter fields is very important, the analysis of the vacuum case is the first thing to address. We 
will come back to this point in the discussion section. 

The paper is structured as follows. In section II we introduce the model together with the equations of motion for 
a FLRW background. In section III, we perform a perturbation analysis for the scalar modes and derive a master 
equation. In section IV, we analyze the master equation and discuss properties of the solution. In section V, we 
consider a toy model of MGM to demonstrate how our findings can actually set bounds on MGM. Section VI is 
devoted to the study of special cases where the structure of the master equation for high k modes differs from the 
general case. In section VII, based on our findings obtained from the study of the scalar perturbations, we discuss 
various implications which must be taken into account for building up a sensible MGM. We also comment on the past 
studies where the scalar perturbations for the MGM were studied. Finally, in section VIII we present our conclusions. 



II. MODIFIED GRAVITY MODELS 



We consider the modified gravity action given by 



-9fiR,G), 



(1) 



where R is the Ricci scalar and G is the so-called Gauss-Bonnet term defined by 

g = i?2-4i?^,ff"' + i?^,„^ii""'"'3. (2) 

Since we want to study the general theory of cosmological perturbations for the modified gravity action ([T]), we do 
not assume any particular functional form of f{R,Q). 

Although ll} is the basic action we consider, for convenience of the actual analysis, we mainly use a different action, 
equivalent to ([1]), which is given by 



S 



Ml 
167r 
Ml 
167r 



fix, a) + iR- X)FiX, a) + ig- a)^{X, a) 
RF{X,a)+g^{X,a) - V{X,cr 



where A and a are auxiliary fields and 



F(A,a) 
V{X,a) 



XF{X,a)+cr^{X,a) 



/(A,a). 



(3) 

(4) 
(5) 



By the following way, we can verify that the action ([3]) is equivalent to ([T]). From the variation of S with respect to 
A and a, we have equations for A and a which are given by 



where Fa = dF/dX, F„ = dF/da and C„ = 
are independent and A and a are given by 



(i? - X)Fx + {g- C7)F, = 0, 
(i? - X)F, + {g- = 0, 

d^/da. If a combination Fa^ct 
A = i?, 

a = g. 



(6) 
(7) 

F^ does not vanish, the two equations 



(8) 
(9) 

Eliminating A and a in the original action by using these results, we find S in Q reduces to ([T]): the equivalence of 
these two actions also holds for the equations of motion. If the combination F\^cr — F^ vanishes, ^ and ([7]) are not 
independent from each other. This does not mean we can not eliminate A and a from the original action. Because ([6]) 
is the condition that the action remains the same under the variation, all the possible pairs of A and a that satisfy 
dS]) give the same action. We can put any (A, a) we want into S as long as (A, a) are the solutions of Obviously, 
X = R and a = g are the solution and we can put them into ^ to find that it recovers the original action P]). 



A. Equations of motion 



The equations of motion for g^^, in the presence also of a matter component with stress-energy tensor T^i,, are 
given by 



where is the effective energy momentum tensor defined by 



The equations of motions for A and a are instead given by A = i? and a — G- 



(10) 



(11) 



B. Background dynamics 

We assume that the background spacetime is a flat FLRW universe whose metric is given by 

ds^ = -df2 + a^(t) dxi dx\ 
where the indices ... are raised by Sij. Then the background components are 

Go Q Tj2 

— ^oti , 



from which we have, in vacuum 



SHF + 12H'^( -\V 
2HF + SH-i + m'^'i + F -\V 



R = 6{2H^ + H), 
g = 24H^{H^ +H), 



3H^ = - 

H = 



i y - iHF - l2H^i 



1 



2F + 8iJ^ 
A = 6{2H^ + H), 
a = 24:H^{H^ + H). 



-F + HF-AH^i^-Hi) 



(12) 

(13) 
(14) 

(15) 
(16) 

(17) 
(18) 

(19) 

(20) 

(21) 
(22) 



III. SCALAR PERTURBATION 



A. Perturbed metric 



We consider scalar perturbations around the metric Eq. p2[) . Wc write the perturbed metric as 



ds' 



-{1 + 2a) dt' - 2a{t)d,f3dtdx' + a^(t)(% + 2(f>6,j + 2d,d.j-f) dx' dx^ . 



For later convenience, we define x by the following equation, 

X = a{P + aj). 



(23) 



(24) 



This represents the shear potential of the unit vector normal to Ei, where is a time-like hypersurface of constant t. 
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B. Gauge transformation 

There are degrees of freedom of choosing Sf. Changmg from St to Et corresponds to a time coordinate transfor- 
mation: t ^ t + T{t,x^). Under this transformation, the perturbation variables transform as 



a^a-t, (25) 

= 0- HT, (26) 

X^X-T, (27) 

SF = SF - FT, (28) 

(5C - - iT. (29) 



C. Perturbation equations in general gauge 



The perturbations of the Einstein tensor can be written in Fourier space as 



(5Gj] = 2 ( 
6G° = 29, 
6G 



Ha), 



3H^ + ^Hx 



y^SG'^ = a-' (^d^d, - i<5;. a) (^ + Hx~cI,- a). 



(30) 
(31) 
(32) 



It can be shown that under the gauge transformation t t + Tit, x^) the Einstein tensor transforms as 

5Gl = 5G°o + &HHT, (33) 

(34) 

(35) 



= (5G" - 2Hd^T, 



SG) - l-S'SGl, = 6G] ~ \5]5G\. 
-'3 ■' 6 ■' 



We do not consider the trace part of because the perturbation equations for the trace part can be derived by the 
combinations of other equations and do not bring new informations. We find that the traceless part of SG^j is gauge 
invariant. 

The perturbations of given in Eq. pT|) can be written as 



<5S 



,5S" = 



1 



h2 p 

3(F + 12i/2^)0 + 8^H^(j) - 6H{F + SH^^a - -^{F + 12H^^)x 



3H{SF + iH^SC) + ( -3{H'^ + H 



(SF + AH^Si) 



—d,[-S'F + HSF - iH^SC + 4:H^6C + (F + 12H^^)a - 8iJ^0] , 
F 



a^F 



1 



^S)Aj |4i/e« + 4^0 -[F + A{H' + H)i + AH^] x 
4H^X + SF + 4.{H^ + H)6^ 



(36) 
(37) 

(38) 



With the help of the background equations, we find that under the gauge transformation t t+T{t, a;*), transform 
as 



(JEg = 6T.°a + 6HHT, 
ST. 



(39) 
(40) 

(41) 



The explicit verification that both G(^ and transform in the same way under the gauge transformation supports 
that the derived perturbation equations are indeed correct. 
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Collecting these results, the perturbation equations, in Fourier space, are given by 



1 



■ ■ P ■ 

3(F + l2H^^)(j> + 8^H^(f> - 6HiF + 8H^^)a 



— {F + 12H^C)x + 3H{dF + AH^SO 



-3{H^ + H) 



{SF + AH^6^) 



(j)- Ha 
X + Hx- 

S\ = -2 
da = -8 



2F 



H5F -5F - AH^5£, + AH^5£, + {F + \2H^i)a - SH^cj) 



1 



AH^X + 6F + 4(ij2 + H)S^ 



1.2 l2 l2 

6{2H^ + H)a + -^{x + 2Hx) - 2^0 -^a-iij)- 12H^ + 3Ha 

12(ij2 + H)H^a + 2^H{H^ + H)x + -^H^X - ^H^(l> + SH^a 
a a 

. . 1.2 , fc2 ' 

- 6H12H^ + HU - 2^(H^ + HU - -^H^a 



(42) 
(43) 

(44) 
(45) 

(46) 



D. Gauge invariant variables and the master equation 

It is convenient to analyze the perturbation equations in terms of gauge invariant variables. Therefore, let us define 
the following gauge independent combinations of fields 



5r = FS£_-^SF, 

H{SF + AH'^S^) 



F + 47?2^ ' 

5F + AH^Sf , 

+ - Y 

2{F + AH^) 2 



FH-F 



F + AHi 



(47) 
(48) 

(49) 



We will see that these three fields are enough to understand the behaviour of the scalar perturbations in the most 
general theory of gravity /(i?, Q). This field definition cannot be applied on maximally symmetric backgrounds, where 
H = Hq, and F = ^ 0. These backgrounds must be studied as a special case, which will be done later. In general, 
one can solve equation for a and substitute it into equation (|^^ . Afterwards one can use the fields introduced 
before, to change (|42|) into 



— ^ Ai{t) <i> + A2{t) ST. 



(50) 



where Ai^2 are functions of the background only and are defined in the appendix. It must be noticed that, in order 
to write down this equation, we made use of the equations of motion in order to replace F in terms of the other 
backgrounds variables. Along the same lines one can rewrite equation (j44p into the following form 



^ = Asit) 6r + A4{t) + Ar,{t) $ . 



(51) 



We need at least another equation in order to make the equations closed. For this aim, we use then the definition of 
the field ST and write the following formula 



ST = {F^, - iF,)5a + [F^x - ^Fx)dX = {Fx - F^) [Xda - aSX) 



(52) 
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It is possible to rewrite both F\ and in terms of F, ^, and Fa- as follows 

(54) 

a 

Therefore we can put equations and (|46|) into equation ([5^ to write another equation for ST and the other fields. 
The equation that one obtains has the following form 

Pi {t)ST+ p2 (t) ^+p3it)ST + — [p4 (i) <5r + * + p5 (t) * + P6 (0 + P7 (i) * = , (55) 

where the pi are all functions of time only. This complicated equation can be simplified as follows. From equation 
(|5T|) one can solve for ^1* = ^'(5', $, ST). By using equation ((50|) one can find $ = $(^1', (5r). By differentiating once 
this last equation and by replacing and F (this last one, by using the background equations) one can also find 
$ = $($, ^P, (5r, ^r). Inserting this last equation into ((55|) to eliminate $ and also replacing $ and ^E* with expressions 
that do not contain time derivatives of the fields, we have an equation given by 

Sr=^[Ae{t)<i> + Aj{t)^i!], (56) 

where ^6,7 are given in the appendix. We find that 6T automatically disappears from the equation, which is crucial 
to derive the closed second order differential equation. As we will see later on, some gravity models identically yield 
i5r = 0. In these models, we can verify that Aq and Ay also identically vanish and (|56[) does not give any information. 

Now, it is possible to find a closed second order differential equation in time for the Fourier fields or ^, or ST. 
For example, let us use equation into equations ([5T|) and ([501) obtaining 

^[l-A2 it)Ar it)] VI/ = Ai (t) $ + A2 (t) Ae (0 ^ 1> , (57) 

* = Asit) Aeit) $ + A^it) Ay{t) ^ + Ai{t) * + A^it) $ . (58) 

Then let us solve ([57| for 'J, finding vp = vl/($^ <|)). Then we can use this result into equation ((55|) to find an equation 
of the form 

^ Ot[a3g(i)<i']+i?i(t)^$ + B2(t)^* = 0, (59) 



a3 Q{t) 



where Q, Bi, and B2 are defined in the appendix. The field <& indeed seems to be a "good" field to study in the sense 
that its equation of motion has a relatively simple dependence in time and k. 

Eq. (j59p is the main result of this paper. The scalar perturbations have only two independent degrees of freedom, 
which are $(io) and $(io)- AH the informations of the perturbation behavior are completely determined by ([55)1 . To 
see this, looking at equation ([57]) and ([55)1 . one can write both \1/ and ST as linear combinations of <I> and <f>. As for 
the metric perturbation variables themselves such as a, x and </>, because of their gauge dependence, they cannot be 
written in terms of 4>, and dV which are gauge invariant. But this does not mean that we need further knowledges 
which are not contained in vp and 6T. It merely means that only the gauge invariant combinations out of {a, Xi 0}j 
or equivalently {a, x, 4>\ in a specific gauge, are expressed by <&, and ST. As an illustration, let us consider a gauge 
where SF = — 4if^(5^, which we call the MGM Gauge (MGMG). There are no remaining gauge degrees of freedom in 
this gauge. In this gauge, we have </) = <!> from (|^5)) . By using the remaining two equations (|T7|) and (|m|) . (5^ and x 
are uniquely expressed in terms of $, vp and ST. After that, a is given by (j44p and finally all the metric perturbation 
variables are determined. Once the metric perturbations are determined in a specific gauge, those in other gauges are 
simply obtained by gauge transformation. 

Therefore, the knowledge of $ and <!> is enough to understand the behavior of the metric perturbations. 

The master equation (j59p contains a term proportional to fc"*, or equivalently a term of fourth order spatial derivative 
in real space. This term does not vanish in generic /(7?, Q) models, nor is a spurious result due to a bad choice of 
fields/gauge: the field <& is gauge invariant, and this fc^ behaviour would still be there by studying either ST or vj/. 
Furthermore in the MGMG $ has a simple meaning which can be directly related to experimental data. Therefore, 
for the general theory, we don't have a standard wave equation and we expect that the perturbations propagate in 
space in non-trivial ways. 
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It should be noted, however, that B2{t) identically vanishes if -Fa Co- ~ = 0. Interestingly, most of the modified 
gravity models concerned by the past Uterature, such as f{R) gravity and R + f{Q) gravity models, belong to these 
special cases. In these special cases, the propagation properties of the perturbations deeply differ from those of the 
generic cases. We defer to study the special cases in section VI. In this section, we consider the generic cases where 

The extreme complexity of Bi (t) and B2 (t) requires numerical calculations to solve exactly the differential equation 
and to see how the perturbations evolve. However, studies for limiting cases still allow us to obtain important 
informations such as instability of the perturbations, which will be done in section IV. 



E. Action of the perturbations 



By expanding the action ([3]) up to second order in the perturbation variables and eliminating all the auxiliary fields 
using the equations of motion, we find that the action in terms of can be written as 



TT 



(60) 



From this action, we correctly recover the master equation (|59p. We find that the sign of the kinetic term in the action 
is equal to that of Q. We will then call that mode for which its kinetic energy has a negative sign, that is Q < 0, a 
ghost. Typically, if a ghost is present and coupled with other normal fields, we can expect that the vacuum decays 
into metric and normal fields, because the ener gy c onservation allows such a process to happen. Therefore, one would 
expect strong bounds from the vacuum decay [TOj. Here we give the condition Q > as a necessary one in order to 
have a theory without ghosts degrees of freedom. However, we are not interested in the details of the quantization 
procedure which is outlined in p^ . Although the quantization mentioned in (69j only treats the case B2 = 0, the fc"* 
term will only change the dispersion relation LUk(t) for each decoupled mode. Therefore, this change will only lead to 
an explicit expression for the solution which is, in general, different from the B2{t) = case. Anyhow, we will 

leave this point in detail for a future project, when matter fields will be introduced. In summary, the first bound we 
will consider is the one coming from not having ghosts in the theory, that is requiring a positive Q. 



IV. STUDY OF THE MASTER EQUATION 



A. Long wavelength Hmit 

Let us first study a case where the wavelength of the mode we consider is much larger than any typical length scale. 
In this case, neglecting terms that are quadratic and quartic in k/a in (j59p yields 

^ at[a=^g(t)$] ==0. (61) 



Q{t) 

We can immediately integrate this differential equation and the result is 

/■* dt' 

where Ck, are independent of t, being integration constants for each k. 

4" can be determined by substituting the above solution into ([57|) . To the leading order, 'J is given by 



1 



Depending on the leading power of k for both Ck and D^, each term on the right hand side can be dominant. 



B. Short wavelength limit 



Let us next study a case where the wavelength of the mode we consider is much smaller than any typical length scale. 
In the large k limit, the time scale for the change of $ is much smaller than that for the change of the background 
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quantities. Therefore, we can use WKB approximation to obtain the solutions. The solutions of ((59|) . under the WKB 
approximation, can be given by 



w c+exp 



+c_ exp 



•,(+) 



2J+^ -i{a^Qy/{a3Q)^ 
,■/,(-) 



24-)-i(a3Q)7(a3Q)^ 



(64) 



where c± are constants and uj^^' are the roots of the following quadratic equation 



fc4 fc2 

oj - B2 - -^Bi = 0. 



In the large k limit, w^^'' are given by 




2 a^Q 



(65) 



(66) 



We see that if B2 is negative, then w^^'' are pure imaginary. Therefore, $ grows exponentially in time. The growth 
rate increases in proportion to k^. The smaller the wavelength of the mode is, the smaller the time scale of the 
instability is. One may guess this instability is due to a bad choice of perturbation variables and there might exist a 
special gauge where all the metric perturbation variables do not exhibit such an instability, i.e. where the deviations 
from the FLRW universe remain small. If this were possible, written in terms of the metric perturbation variables 
in that gauge, should not show an exponential growth, which is inconsistent with the exponential growth of <f> we 
have just found. Therefore this instability is not a due to a bad choice for the perturbation variables and the FLRW 
universe is indeed unstable on small scales. Due to this instability, the perturbations will grow until they become 
C(l) where the linear perturbation theory no longer works. We do not know what happens if the system goes into 
the non-linear regime and will not consider it furthermore. 



If B2 is positive, then the leading term of Luj^ 
velocity is given by 



(±) 



is real and the perturbations propagate. In this case, the group 



Vg{k) 




(67) 



At any time to, this velocity exceeds the speed of light for those modes above a critical wavenumber kc which is given 
by kc = a{tQ)/[2^ B2{to)]. The propagation of the short wavelength modes inevitably becomes superluminal up to 



the cutoff of the theory w„ 



/_B2Acutoff- The smaller the wavelength of the mode is, the larger the propagation 



speed. For a given fc, the speed of propagation is background dependent. Because of this, if Vg reduces quickly with 
time, for some backgrounds, it may be possible that the friction term (the third one in the right hand side of Eq. I66p 
gives the strongest contribution to the stability of each mode, at least for t — > 00. It still remains a lack of consensus 
among researchers if the propagation speed larger than the speed of light causes crucial problems that immediately 
make us give up considering such models as a possible modification for a consistent theory of gravity. More details 
are given in the discussion section. 

We will also impose Bi to be positive as there is always some intermediate range of fc's, for which the Bik^ term 
will be dominating over the B2k'^ one. 

To summarize, for the general f{R, G) MGM models, except for those special cases where Fa^o- — = 0, the small 
wavelength modes inevitably either suffer from strong instability or acquire superluminal propagation. 



V. APPLICATION TO A TOY MODEL 



Just for the purpose of illustrating how the general formalism we have developed in the previous section can be 
applied to an explicit MGM, let us consider the following toy model, introduced in ([l]), whose f{R,Q) is given by, 

fiR,g) = R+ (68) 
R'' + ay 
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where M is a constant of mass dimension and a is a dimensionless constant. In this model, we have that 



dwWydim) ^'{m + agf 



(69) 



which docs not vanish in general. Therefore, this model docs not belong to the special class of MGM, and its master 
equation will have the fc^-term. 

It is known that the background equations can have power-law solution such as a{t) oc for four different values 
of p. They are given by 

PI = '^tr^f^"^ (70) 



4(6 + a) 
51 + 10a -;/(«) 

P2 = . s , (71) 

4(6 + a) 



6 + a + \/6a + c? 

= — WT^ — ' ^ ^ 



6 + a — \/6a + c? 
2(6 + a) 



where y = y{a) = \/1521 + 840a + lOOa^. If the radicand becomes negative, which indeed happens for some values 
of a, the corresponding power-law solution will not exist. 

From now on, we only consider the case p = pi. Analysis for other pi can be done in exactly the same way. From 
the background equation, we find that the power-law solution is an attractor if the following inequality 

273 + 50a + 3y(a) 

6+a ' ^ ^ 

is satisfied. If not, the universe cannot settle down to this power-law solution. From the definition of y{a), we find 
that y{a) is real for a < -3(14 -h 3V3)/10 « -5.76 and for a > 3(-14 + 3V3)/10 « -2.64. For such ranges of a, the 
power-law solution becomes an attractor for a < —507/80 « —6.34 and for a > —2.64. Moreover, such an attractor 
is an accelerating one for a > —2.64, whereas it is a decelerating one for a < —6.34. 

Using the formulae given in the appendix, the expressions for B2 and Q at late time (t — > 00) are given by 

^ _ 135(26364 - 391y) + 16a{60(3660 - 49?/) -I- a[70716 - 631?/ -t- 10a(907 + 40a - 4?/)]} 
^ ^ 45(15 + 8a)(57 + 8a)?/ ' ^ ' 

jj2j^ ^ -16a^(6 + a)(51 + 10a + ?/) 

^ 3[8a(15-f 2a) -f 3(39 4- y)]{351(39 + y) + 2a[5301 + 87?/+ 10a(129+ 10a + ?/)]}' ^ ' 

^6q ^ 10125(15 + 8a)(57 + 8a)?/ 
4s(a) 

s(a) = -4108137345(39 + y) - 2a{-9477(-18610644 + 44279?/) + 10a [-297(-46531983 + 545000?/) 

+ 320a (27(579515 - 6608y) + 5a[581040 - 4563?/ + 10a(5403 + 200a - 20?/)])] }. (78) 

Here we have multiplied B2 and Q by and to make them dimensionless. Graphs of these quantities as functions 
of a are shown in Fig. [TJ 

Here we find that the region of a where all of the quantities become positive is —2.64 < a < —2.25. This range 
is outside the allowed range coming from different constraints regarding the fit to Supernovae Type la data and the 
viability of the cosmological evolution from radiation domination up to today [g^, [7l|. Therefore, according to our 
results, for these models, in the a-interval allowed by theory/experiment, the background attractor pi will be in 
general unstable. In fact, outside this region, at least one of {Bi,B2,Q} becomes negative. The asymptotic value of 
H^B2 for a —00 is —1/9. At first glance, this is incompatible because the action ((68)l in the limit a ±00 reduces 
to a form f{R,Q) R + M^/{aG) which belong to the special case. The origin of this gap comes from the naive 
expectation that in the denominator of ([55)) will be much smaller than aQ for |a| ^ 1 and can be neglected. This 
expectation is true for positively large a. For negatively large a, pi becomes 9/(4a) and R^ in the denominator of 
(j68p remains of the same order of magnitude as aQ. Therefore, as long as this this background solution is concerned, 
the perturbation behaviors in the limit —a 3> 1 do not reduce to the special case. 
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FIG. 1: Plot of Q,Bi,B2 for one attractor of the toy model considered. Bi becomes positive for three regions of a: (i) 
a < -8.19, (ii)-2.64 < a < -2.25, (iii) a > -1.875. H'^B2 becomes positive only for -2.64 < a < -1.875. H^Q becomes 
positive only for —2.64 < a < —1.875. 



VI. SPECIAL CASES 

In the previous section, we found that the perturbations propagate in non-trivial ways on small scales due to the fc^- 
term in the master equation ([5^ . However, there are special cases where B2, the coefficient in front of k'^^, vanishes 
identically. In these cases, we do not have the exotic behavior such as a strong instability or a superluminal propagation 
observed in the generic models due to the fc^-term. Instead of the fc^-term, the fc^-term becomes important. In this 
section, we study these special cases. 

From the explicit expression of B2 given in the appendix, wc sec that the special case happens if either TJ = or 
F\^^ — = 0. Therefore, dc Sitter space for any f{R, G) MGM belongs to this case. Wc can easily verify that both 
R + f{G) and f{R) models for any background solutions also belong to this case. As wc will explain later, there are 
other infinite numbers of f{R, Q) models that, irrespective of the background solutions, yield B2 = 0. We will provide 
a systematic method to find those special gravity models and will list some explicit forms of /(i?, Q). 

Let us comment that vanishing B2 is equivalent to 5V = 0. To see this, let us first assume B2 = 0. This implies 
that 7? = or F\£^„ — F^ = 0. If H = 0, then A and a are also zero. From ([5^ . we immediately find that ST is 
also zero. If F\^cr — F^ ^ 0, (|52p tells us again that ST = 0. Therefore, B2 ~ implies ST = 0. Next let us assume 
that (5r = 0. This implies that either of Fx^cr — F^ = or X 6a — a 6X = 0. In the first case, it is obvious from the 
expression of B2 given in the appendix that B2 vanishes. In the second case, both A and & must vanish, which implies 
H ^ 0. This again gives B2 = 0. Therefore, ^2 = is equivalent to JF = 0. 

Since there are two independent conditions for the models to belong to this special case, we will study them 
separately. 

A. Special cases 1: H = 

If the background space-time is de Sitter, we have A = ct = 0. Since F and ^ are both functions of A and a only, 
we also have that F = £^ = 0. Except for the scale factor, which varies like a{t) oc e^*, all the other background 
quantities appearing in the perturbed equations arc constants. Keeping this in mind, the perturbed equations (j42p . 
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(|43)) . dHD, (gSl) and ^ become 



2F 



3H{6'F + 4,H'^S0 



{SF + AH^SO 



^- Ha = — [HSF -6'F- AH^S^ + AH^S^] , 
X + Hx-<t>-o^ = ^[5F + AHH(\ , 



5\ 
da 



12H^a + — (x + 2Hx) - - —a - 30 - 12Hcj) + 3Ha 

l2H^a + 2—H^x + -^H^X - SH^(f> + ^H^a - 12H^cj) - 2—H'^cj) - —H'^a 



(79) 

(80) 
(81) 
(82) 

(83) 



These equations can be rewritten in terms of gauge independent fields. We can adopt the same field 4" as in the 
general case (but evaluated on the de Sitter background), however $ is not well defined for a de Sitter background. 
We will choose then another field $ds defined as 



SF + AH^S^ 
2F 



(84) 



This combination is gauge invariant on de Sitter and reduces to </> in the Newtonian gauge (x = 0). Therefore the 
equations of motion can be rewritten as 



* = 0, 

$DS = -3-ff$DS 



'i-AH' 



F 



3{Fx + SH^F^ + 16m^„ 



DS ■ 



(85) 
(86) 



where SF and S^ are not independent from each other as, on this background, we have Sa = AH^SX. Eq. ^ is equal 
to the Klein-Gordon equation in de Sitter space-time with an effective mass 



-AH-" 



3(Fa + 8H^F^ + l&H^ia) ■ 



(87) 



If JTi^ff is negative, then the long wavelength modes are unstable (tachyonic instability). For any f{R,Q) gravity 
models, the short wavelength modes arc stable and the sound velocity is unity. The last quantity we need to find for 
this special case, is the expression for Q for ^us- Its action can be written as 



Mi 



dtdPxQ 



^ *DS - ^ ^ (V$DS)' - \ mis *DS 



where 



(89) 



Therefore F must be positive in order to avoid ghost-like degrees of freedom to propagate. 



B. Special cases II: FaC<7 - = 

Before wc study the perturbation equation, let us try to find what kinds of models belong to this case. 
Because F and ^ are functions of A and a, we have 

SF = FxSX + FM, (90) 
(5e = FJX + S^Ja- (91) 

Here we have used the identity ^\ ~ F„ in the equation for d^. Because i^A^cr — F^ is the determinant of the matrix 
of the linear transformation above, the condition Fx^cr — F^ ^ tells us that SF and S(, arc not independent. This 
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implies that F is a function of ^. In other words, F and ^ are related to each other by a single variable ip like 
F = F{ip), ^ = S,{ip). Such modified gravity models are equivalent to the action 

/\/2 r . . _ 

^ - Ye? y [^(^)^ + ^(^)^ - (92) 

where F, ^ and [/ are arbitrary functions. To see that this action is equivalent to the f{R,Q) action for the special 
case, let us first take the variation of (j92p with respect to if. After the variation, we obtain an equation given by 

dp df dU 

dip dip dip 

By solving this algebraic equation with respect to ip^ we can write as a function of R and Q, i.e. ip — p{R,Q). 
Substituting this solution into ([5^ . we find that it reduces to the f{R, G) action. Then we can calculate F and ^ as 

F^^=F{p^iR,g)), (94) 
( = ^=i{'PiR,G)). (95) 



Therefore, F and ^ are equal to F and ^. From these equations, we can easily verify that this model belongs to the 
special case, since 

dRdg^dGdR^ ^ ' 

In this way, we can construct f{R, G) models that belongs to the special case. Though finding analytic solution 
of (|93p is difficult in general, we can in principle find as many models as we want by doing the above procedures for 
various F^ ^ and U. In this paper, we provide only a few simple f{R, G) models by using the above method. The first 
example is a case where ^{(p) = c (c is a dimensionless constant). In this case, ([M)) tells us that (p depends only on 
R. Therefore, the corresponding modified gravity model is f{R, G) = f{R) + cG- Note that the term cG only adds a 
total derivative in the action and does not contribute to the equation of motion. Therefore, this model is equivalent 
to the f{R) gravity model. The second example is the opposite of the first case, i.e F{ip) = c. The corresponding 
modified gravity model is f{R, G) = cR + f{G)- The third example is a case where F{ip) = A'P(^{ip) (M is a constant 
of mass dimension). In this case, ([M)) tells us that ip depends only on the combination R + Q/M^. The corresponding 
modified gravity model is f{R,Q) = f{R + Q/Al"^). We list these models in Table H] For the special models, the 
perturbation equation ([59|) can now be written as 

^ ^t[a''Q{t)^+B^{t)%^^Q, (97) 



a3Q(t) ' ""a 

where Bi in these cases is significantly reduced to a much simpler expression than the general case and is given by 

^ ^ 16H^[F + 4(g2 + H)i] + F[3F + 4(3g^ + AH)i] - 4{F + AH^)^ 

'i{F + AHi){F + AH^i) 

In particular, for f{R) gravity models, it reduces to 

Bi{t) = l. (99) 

Therefore, we correctly recover the well-known result that the propagation speed is equal to the velocity of light in 
f{R) models. 

For cR + f{G) gravity models, the propagation speed was calculated in [6^. In this case, according to our formula 
((98)) . Bi{t) reduces to 

^iW = i + ii- (100) 

If the universe, in vacuum, accelerates more slowly than a{t) cx i^, Bi becomes negative and the FLRW universe is 
unstable on small scales. Therefore for these theories, acceleration (faster than is a condition for vacuum stability. 
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Special f{R,G) models 


remark 


f{R) + g/M^ 


£,{ip) = MiConstant of mass dimension 


cR + f{g) 


F{ip) = c, c:Dimensionless constant 


f{R + g/M^) 





TABLE I: We list some of the special f{R,g) models for which B2{t) identically vanishes. In the first case, different values of 
M give the same modified gravity model because a linear term in g in the action is a total derivative and does not contribute 
to the equations of motion. 



On the other hand, if the universe undergoes super-acceleration H > 0, the propagation of the perturbations becomes 
superluminal. 

For f{R + G/M'^) gravity models, Bi{t) reduces to 

^iW = i + ^^^^. (101) 

Therefore, similarly to the cR + f{G) models, the propagation becomes superluminal if the universe undergoes super- 
acceleration. 

VII. DISCUSSION 

So far we have considered the scalar modes. In the literature [6^, [6^, it has been claimed that the cosmological 
scalar perturbations for a general Lagrangian of the kind R) + £,{(/>) G — Ui^cj)), were studied^. However, we think 
that the scalar perturbations of this Lagrangian were not studied there. In fact, the authors considered only the 
following Lagrangian R + £,{(j))Q — U{(j)), which actually reduces to the special case discussed above. By doing 
so, in the literature, the fc^'-term has always been neglected, and the most general MGM was never then fully studied. 
In fact, the novelty of our analysis resides in the study of the general case where the modes do get a non-trivial 
modification for the dispersion relation. 

As for the vector modes, they are not important since, as already stated in [s^, they do not propagate. We have 
explicitly verified that this remains true even for the general cases for which Fx^o- — ¥^ ^ ^'^'^ H 0. The tensor 
modes instead, do propagate, but, differently from the scalars, we checked that there is no difference between general 
and special cases, and that their evolution equation coincides with the one given in [69j . More in detail their equations 
are 

1 - 
^T;— dtia^'QTTC^,,) - W^Cf,. = 0, (102) 

where 



iTT 



2 



F + 2Hi, (103) 
c|j, = . (104) 

However, they still give an important contribution in order to set constraints for these theories. To summarize, we 
can then set the following points in order to understand the vacuum structures for the MGM. 

• In the general case one must set the following conditions for the tensor modes Qtt > 0, c^rp > 0. The scalar 
modes instead have a more complicated structure as they will acquire a non-trivial dispersion relation. We 
found, at least for the long and short wavelength modes, that Bi > and B2 > are the necessary conditions 
to avoid instability. Furthermore, in order to avoid ghosts, one needs to set the condition Q > 0. 



^ Eliminating the auxiliary field ip by using the equation of motion, we find that this model is equivalent to the general f{R, Q) gravity 
models. 
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• The backgrounds for the general MGM, which are not de Sitter, if all the previous conditions are satisfied, at least 
in the scalar sector, will always have modes with superluminal propagations. In fact, for the short-wavelength 
ones, independently of the background, their propagation speed grows linearly with their wave number fc/a, up 
to the cutoff of the theory. The presence of superluminal modes might be present also for the tensor modes, 
but this is only a background dependent feature, i.e. the propagation speed is independent of k. Recently, in 
the literature, superluminal modes in cosmology have been discussed in the context of different theories, such as 
k-essence models. Whether or not these modes represent a real issue on a background which explicitly brakes 
Lorentz invariance is still matter of discussion and a clear and definite opinion shared by everybody on this 
issue is still not achieved [zl, [zl, [zl, [ll, S S H Izl, S HH, HI . Some problems may arise from having a 
superluminal modes interacting with matter fields, either indirectly (by looking for example at the evolution of 
the clustering Sp/p) or directly (such as the production of a kind of Cerenkov radiation), but this would imply 
to go beyond linear perturbation theory. 

• The behaviour of the scalar modes can be used in order to distinguish among different MGM. In fact, we have 
shown that the modified gravity models which possesses second order differential equations can be divided into 
two categories 

1. General Gravity Models, for which — Fx(,cr 7^ 0. The scalar modes have two degrees of freedom, but 
their equation of motion is non-trivial in the sense that it has also a V^4> term. This terms, at least for 
the short-wavelength modes, changes the dispersion relation such that the group velocity of these modes 
becomes proportional to k, i.e Vg = 2^/B2k/a (if B2 > 0). Therefore the theory, unless the vacuum is 
de Sitter, will always have superluminal modes. This is a new feature of these model which, up to our 
knowledge, was not considered before in the literature. It should be noticed that this feature is not a 
spurious one, in the sense that it can be gauged away. In fact, the analysis has been done through gauge 
invariant fields. Besides, no matter which gauge invariant field is used to study the large k behaviour, 
exactly the same B2 coefficient appears, as the parameter in front of the fc'*-term. 

2. Special Gravity Models, for which — -Fa^o- = 0. The f{R) and R + f{Q) theories belong to this class 
(however general fi{R) + f2{G) does not), but they are not the only ones. In fact we have found that there 
is a larger group of theories which all share the same feature B2 = 0, which distinguishes them from the 
General case. This means that the equation of motion for the modes is similar to a standard wave equation, 
where speed of propagation is background dependent (except for the ,f{R) theories for which = 1). As 
far as we know, these are the models mostly studied in the literature, but actually they represent a special 
case of the whole class of MGM. 

3. The de Sitter background is a "good" background for all modified gravity theories, provided that it is 
stable, that is if the scalar field has an effective positive squared mass, and provided that F > 0, in order 
to remove ghosts degrees of freedom. 

In this paper we have only considered the vacuum case. Although the presence of matter is important, the vacuum 
itself has presented a rich structure due to the presence of the fc^-term. Thinking on how matter could change the 
whole picture, it is difficult to imagine that this fc^-term discontinuously disappears if we add standard matter fields 
(which, in the action, are expected to contribute only to the fc^-terms) into the theory, because this k^ dispersion 
relation is solely due to the modification of gravity and not due to the nature of matter. In this sense, the essence of 
the k'^ propagation, already present in vacuum case, clearly supports the importance of studying the vacuum. 

The next obvious and important extension of this work is to add matter fields into these MGM, that we will study 
in another project of ours. Although conclusive results will be given in our forthcoming paper, we can naively expect 
the following effects on the propagation of modes due to the existence of matter. First, the master equation for 
the scalar perturbation of gravity has time dependent coefficients determined by the background dynamics. Since 
the background dynamics changes if matter fields are present, one obvious effect, by including the matter, is that 
those time dependent coefficients in the master equation will change. Second, we are adding new degrees of freedom. 
Therefore, we expect that one variable is not enough to specify the perturbation behavior. We will obtain at least 
two coupled evolution equations for the perturbations. 

VIII. CONCLUSIONS 

We have studied the structure of the scalar cosmological perturbation for the most general classes of MGM, /(i?, G) 
which do not have spurious spin-2 degrees of freedom. Indeed we have discovered, for the general case, a non-trivial 
propagation for the modes that has not been studied before. In fact, in the past literature [661. [69j. claims were made 
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about studying the general class of MGM. However, we believe that this study was never actually performed, as the 
authors only discussed a special subcase of the general MGM. 

This new propagation affects primarily the short-wavelength mod es on n on-de Sitter backgrounds, for which their 
group velocity becomes proportional to their wave number, Vg = 2-^/ B2{t)k/a, where B2 is a background dependent 
quantity. This amounts to have a non-standard dispersion relation. In turn, this implies that these models will in 
general have superluminal modes. Besides we have also given the conditions that the same modes may be stable, i.e. 
they do not acquire a negative squared speed, or they do not become ghosts degrees of freedom. Although there is 
still discussion about the presence in cosmology of superluminal modes, we have here a new necessary constraint to 
impose for these MGM, that is < i?2 < -l-oo. 

However, not all MGM share this feature. Indeed there is a subclass (to which the common f{R) and f{Q) theories 
belong) which does not have any longer these modified dispersion relation. 

We strongly believe, that in the process to find (through a phenomenological bottom-up approach) the ultimate 
low-energy effective theory of gravity (and possibly to find at last that it coincides with General Relativity, albeit a 
cosmological constant), this work can be considered fundamental in order to set constraints to the cosmology of these 
models. 

The search for the gravitational action requires indeed the study of local gravity constraints, the study of cosmo- 
logical solutions, but, also in the light of new experiments, more and more the analysis of cosmological perturbations. 
It is necessary to understand the behaviour of these same models in the presence of matter fields, but we leave this 
issue for future work. 
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APPENDIX A: EXPRESSIONS OF Ai{t), (i = 1,...,7), Bi,B2 AND Q 

We give explicit forms for Ai{t). 

A, it) ^ ^(F^AHH? ^ 

2(F + mS) [F + 2H{F + 6H£,)] 

Al{t) = : : ^ (A2) 

{F ^ mi) {F + 2H{F + 6Hi)] 
Asit)^ . (A3) 

^ ^ 2{(F + AHH)[F + 4(g + H^)i] + FH{F + 8HH) + 4g|[j- + H{F + SffQ] + F^H^} 

{F + AHi)[F + 2H{F + 6Hi)] 
^ 16Hi[F + 4{H + H^)i]{F + iH^i) + F[3F + 4{4H + 3H'^)i]{F + 4H^)~4i{F + m^ 

A^tt) = : : : ^ (A5| 

2{F + AHi)^ [F + 2H{F + 6Hi)] 
^ ^ 2H{F + AHH)^6F^ [HjF + AHH) + ^HH[F + 2{H + 2H^)i]] - Pj} 

3H[F + 2H{F + 6Hi)] [2H{FH + 2H^F - ASF^Hm + 8iJ4^) + HH{F + AH^i)] ' 
^ 2H{32FH^i + iF^H^ + 64:H^i^){6F^ [H{F + iH^i) + 4:HH[F + 2{H + 2H'^)i]] - F^} 

^ 3H[F + 2H{F + 6Hi)][2H{FH + 2mF - 48F^HH^ + 8H'^i) + HH{F + 4:H^i)] 

A useful relation among these quantities is the following A2AQ — AiA^At, moreover both Aq and Aj are proportional 
to F^ — F\icr, so that they identically vanish for the special cases. We also have 

AAF + ^Hi)^ (A8) 



2(^2^7 ~l)[F + 2H{F + 6Hi)] 
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Also we have defined 

-A6A2-Al{A4AeA7-ArAe + AeA7) + A2{A4Ae-Ae) + A5{A2A7-l)^ 
^ AAA2A7-I) ' 

' AlAj 3[Fx + 8H^F, + 2H^i,)]{F + 2iH[H{F,+AH^^„)+mH[F, + 2^,{H + 2H^)]]}' 
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